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Abstract 

Field-angle dependent heat capacity of the non-magnetic borocarbide superconductor YM2B2C 
reveals a clear fourfold oscillation, the first observation of its kind. The observed angular variations 
were analyzed as a function of magnetic field angle, field intensity and temperature to provide its 
origin. The quantitative agreement between experiment and theory strongly suggests that we are 
directly observing nodal quasiparticles generated along < 100 > by the Doppler effect. The results 
demonstrate that field-angle heat capacity can be a powerful tool in probing the momentum- 
space gap structure in unconventional superconductors such as high T c cuprates, heavy fermion 
superconductors, etc. 
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Most superconductors behave "conventionally;" electronic excitations are suppressed by 
the BCS energy gap causing the electronic heat capacity, for example, to be exponentially 
small at temperatures well below the superconducting transition. Recently, a number of 
superconductors, cuprates and heavy-fermion metals among them, have been found to be 
"unconventional" in that they exhibit gap- zero (or nodal) points or lines in momentum space. 
Electronic excitations-nodal quasiparticles (nqp)-are then observed at low temperatures, 
giving rise to power-law rather than exponential, behavior. Unlike gapless superconductivity, 
which can occur in conventional superconductors, the Fermi momenta of these quasiparticles 
are restricted to nodal regions of the Fermi surface, giving a strong directional dependence to 
various physical properties. While pioneering work by Salamon and coworkers Q, Q| and by 
subsequent experimenters jHI^I demonstrated the directionality of thermally excited nqp's 
via the thermal conductivity, direct detection has proven elusive. 

The rare earth borocarbide superconductors RNi 2 B 2 C (R=Y, Lu, Tm, Er, Ho, and Dy) 
lad beengenerally thought to exhibit isotropic s-wave pairing via electron-phonon coupling 
, but there is growing evidence that the gap function is hig hly anisotropic or 



has a nodal structure on the momentum surface jj, llfli fllL 12, Q, 14 1. In specific heat 
measurements, a power law behavior was observed in the temperature dependence and the 
electronic coefficient 7 (if), to the extent it can be extracted, follows a square-root field 

nn 

dependence [10, 15]. Low temperature thermal conductivity measurements showed low en- 
ergy excitations that also indicate unconventional superconductivity 0]. Recently, Izawa 
and coworkers reported a field-directional-angle dependence of the thermal conductivity in 
YNi 2 B 2 C that suggested point nodes along < 100 > directions [4]. However, their inter- 
pretation is not definitive because the highly anisotropic nature of the Fermi surface and 
nonlocality could lead to anisotropic macroscopic phenomena even within the ab-plane [lf| . 
Angular specific heat measurements, however, are insensitive to the Fermi surface topology 
and reflect the gap anisotropy through the Doppler energy shift Q). We report here the 
first observation of Doppler-induced angular dependence of the electronic heat capacity in 
the non-magnetic borocarbide superconductor YNi 2 B 2 C. 

In order to study the elusive angular variation in the low temperature specific heat r], Q| > 
we built a probe specifically designed to look at the DOS variation with field angle. The 
probe was designed to fit into a Quantum Design Physical Properties Measurement System 
with a 7 T transverse magnet. The heat capacity of YNi 2 B 2 C was measured by ac calorimetry 
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in which the sample was heated by an electronically modulated laser pulse and the induced 
ac temperature was measured by chromel-gold/iron (0.07%) thermocouple. In order to 
account for the magnetic field dependence of the thermocouple wire, the ac signal from a 
non-magnetic dummy sample was measured vs magnetic field in the temperature range of 
interest. As the heat capacity is field independent, any change with field arises from the 
thermocouple wire and was used to correct the data. The details of the ac technique can be 
found elsewhere 1201. 



ligh temperature flux method using Ni 2 B 
2ll . |22| . The inset of the Fig. 1 shows the 



The single crystal YNi 2 B 2 C was grown by the 
as a solvent. The details are described elsewhere 
superconducting transition at T. The sharp thermodynamic transition (AT/T C < 0.02) 
indicates that the sample is homogeneous. 

The main graph of Fig. 1 shows the field dependence of the heat capacity at 2.5 K. The 
sample was zero-field cooled to 2.5 K and a magnetic field was applied along the (100) 
direction. The open circles describe the heat capacity with increasing field and the crosses, 
with decreasing field. There is a deviation between the two below 1 T, which indicates that 
flux pinning is important at low magnetic fields. The system enters into the normal state 
above 5.2 T along this direction and the lower critical field is less than 0.1 T. The heat 
capacity with increasing magnetic field was fitted by C + b (H - Ho) 1 / 2 (dashed line). 
The excellent square-root fit up to the upper critical field indicates that YNi 2 B 2 C is a 
superconductor with an unconventional pairing symmetry. In an intermediate field range, 
H c i <C H <C H c2 , the square root field dependence comes mostly from nqp's induced by the 
Doppler effect; the core-state contribution becomes important as the field approaches H c2 
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The upper panel of Fig. 2 shows the field-directional angular dependence of the total heat 
capacity at 2 K in 1 T. The field angle a was measured with respect to the a-axis. The 
transverse magnetic field was rotated within the basal plane of YNi 2 B 2 C by increments of 
3° by a computer controlled stepping motor. Repositioning the sample by 33° relative to 
the thermocouples and apparatus shifts the 4-fold pattern, but not the 2-fold contribution, 
by that angle. The total heat capacity consists of constant, 2-fold, and 4- fold contributions: 
Ctotai{oi) = Cq + C 2 (a) + C 4 (a). The field- independent constant C is due to nonmagnetic 
contributions such as the lattice heat capacity and thermally excited nqp's, and is determined 
experimentally in C vs H data. The 2- fold contribution C 2 (a;) comes from our experimental 



3 



setup and has a functional form of C2Cos2a. The Au/Fe thermocouple wire is a major 
source of this contribution but misalignment of the basal plane of the sample against the 
field direction could also lead to the 2-fold component because of the anisotropy between 
in-field ab-plane and c-axis heat capacities. The dashed line in the upper panel of the Fig. 2 
is C2{a)', the 2-fold signal is about 48 % of the 4-fold component at 1 T and increases with 
magnetic field up to 55 % at 4 T. The circles in the lower panel of the Fig. 2 describe the 
4- fold part 6*4(0;), which clearly reveals the angular variation. 

In the mixed state of an unconventional superconductor, the supercurrent flows around a 
vortex lead to a Doppler energy shift, 5u ~ v s • kp, where v s is the velocity of the superfluid 
and is the Fermi momentum of nodal quasiparticles. When the field direction is normal to 
the plane containing nodes, the DOS is the average over the whole Fermi surface, leading to a 
square- root field dependence ^5|. When the field is in the nodal plane, however, the Doppler 
shift has a field-direction dependence as well, 5u ~ ^ sin (3 sin(0 — a) [17[. Here <fi is an 



azimuthal angle of the gap node and (3 is a vortex current winding angle. Eh = ^-^^JitH/ $ 
is the energy scale associated with the Doppler effect. The geometrical constant a is order 
of unity, vf is the Fermi velocity, and $0 is the flux quantum. The variable p = r/R, where 
r is the distance from the vortex core and 2R is the intervortex distance. Working in the 
2D limit, Vekhter et.al. calculated the DOS , N « (Ni + A^)/2, for an in-plane magnetic 
field when uj, <C A and for four nodes at angles a n from orthogonal axes in a plane : 

_ [ + {x = u/E i >l) 

N o \ + 2x2 ) arcsinx + 3xy/l - x 2 } (x<l), 

where i = 1,2, E\ = Eh\ sin(a n — a)\ and E2 = Eh\ cos(a n — a)|. No is the normal-state 
density of states. The directional dependence of the Doppler effect, therefore, leads to the 
oscillation of the DOS where it has minima along nodes and maxima along antinodes. At 
T = K, the oscillation is cusped with a contrast of l/y/2 between the minima and the 
maxima and is independent of the magnetic field. At a finite temperature, however, the 
sharp contrast is washed out and the DOS oscillation now depends on the magnetic field. 

As seen in Fig. 2, our data are well described by the cusped function | sin 2a \ which will 
be used in the analysis of the 4-fold variation. To make a quantitative analysis, we fit our 
data to Ci(h, a) = 04(1 + T| sin2a|). The coefficient C4 and the angular contrast T were 
treated as field dependent fitting parameters. The solid line in the lower panel of Fig. 2 
shows the fit with T = 0.05 in 1 T. The sharp minima along < 100 > indicate that there 



exist gap minima or nodal structures along those directions. That the nodal positions are 
at < 100 > is consistent with the angular thermal conductivity measurement |4j. 

The contrast T depends on the gap geometry and dimensionality of the superconductor. 
The Vekhter et al. model above assumes the Fermi momentum of nqp's to be restricted 
to the nodal plane. Recent band calculations show that YNi 2 B 2 C has a 3D electronic 

n 

structure in spite of its layered structure p. A more realistic calculation, therefore, is 
to allow nqp Fermi momentum out of the nodal planes, which decreases the amplitude of 



the angular variation. Following Won et al. [24j, we use a modulated cylindrical Fermi 
surface and account for the 3D effect by replacing the Doppler related energies Ei and E 2 
by Ef D = E h [sin 2 (a n - a) + cos 2 6} 1/2 and E% D = E h [cos 2 (a n - a) + cos 2 6] 1/2 respectively 
and integrate the DOS over polar angle 9: 

1 f^l, 



N(w,h,a) = — -[N 1 (w,h,a,e) + N 2 (w,h,a,9)]de. (2) 
2tt J 2 

The heat capacity ratio, C(2K,h,a)/C(2K,h,0), calculated from the DOS, is displayed 
as a 3D surface vs magnetic field- angle a (x-axis) and H (y-axis) in Fig. 3(a). In 
Fig. 3(b), we show the experimentally obtained data sets after subtraction of the back- 
ground contributions. Since the in-plane oscillation is due only to the Doppler effect, 
Ci(2K, h, a)/C^(2K, h, 0) corresponds to the theory. We used 0.6, 1, 1.5, and 2 T data sets 
to construct the contour both in numerical and experimental plots. The Doppler energy 
scale Eh was adjusted to get a best fit, giving the Fermi velocity of Vp = ^/v a v c = 1.3 x 10 7 
cm/s which is comparable to other reported values 25]. We used the geometrical constant 
a = 1 and A = 1.76/cbT c . Note that the angle contrast between experiment and theory are 
in good agreement suggesting that we are directly observing Doppler-shift generated nqp's. 
The periodicity of the peaks is it/2 and the oscillations grow with increasing field because 
of the increased Doppler energy. 

Fig. 4 compares the contrast between node and antinode, T(2K,h) = 
[C(2K,h,n/4:) - C(2K,h,0)]/C(2K,h,0) from experiment (open circles) with nu- 
merical values (dashed line). Above h = 0.5, the data deviates from the theoretical 
prediction. One source of deviation is the Zeeman energy splitting fiH, where /i is Bohr 
magneton, discussed recently by Whelan and Carbotte j26j]. At a critical field H c , Zeeman 
and Doppler energy scales become comparable and the DOS variation with field angle 
vanishes. It is unlikely, however, that the deviation comes from the effect because the 
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critical field is ~ 10 2 T. A more likely possibility is that the deviation is evidence that the 
core-state contribution to the electronic specific heat becomes significant in the high field 
regime. When we include the core states, the parameter C4 consists of the core and the 
nodal contributions : C±(a) = C core + C no d e {oi). Since the oscillation anisotropy T measures 
the amplitude of the 4-fold pattern against the nodal part only, we need to subtract the 
core part from the C4, resulting in an increased anisotropy in our data. In order for a 
quantitative analysis, it is essential to distinguish the contributions from the core states 
and the extended states. 

The inset of the Fig. 4 shows the contrast T vs temperature in 1 T. The basal plane 
of the sample was misaligned by 26° from the field direction, which produced a larger 2- 
fold contribution as was expected due to the heat capacity anisotropy between ab-plane 
and c-axis. The temperature dependence of the anisotropy was compared with the 3D nqp 
model that was used in the above analysis. The dashed line is the estimated anisotropy with 
the same Doppler energy scale with vf = 1-3 x 10 7 cm/sec. It explains the temperature 
dependence of the oscillation amplitude very well, corroborating that the nodal quasiparticles 
are Doppler shifted. 

Before concluding this report, we want to discuss the possible pairing symmetries of 
YNi2B2C. Observation of the four-fold angular variation in our heat capacity measurement 
limits the candidates to those with gap minima or nodes with singlet pairing states such as 
anisotropic s-wave, d-wave or mixtures of two or more order parameters. Recently, Izawa et 
al. reported a clear fourfold pattern in the c-axis thermal conductivity at 0.27 K, which sets 
the gap-anisotropy ratio to be much higher than 100 jj]. One can still argue for gap minima, 
but a more natural explanation will be the presence of nodes in momentum space. We note 
that the field-angle heat capacity will have 4-fold pattern for both d-wave and anisotropic 
s-wave gap with nodes. Considering the significance of the pairing symmetry to understand 
the superconductivity in YM2B2C and further to understand the interplay of magnetism 
and superconductivity in RNi2B2C (R = Tm, Er, Ho, and Dy), other measurements such 
as phase sensitive 7r junction SQUID or tunneling measurements are required to distinguish 
those anisotropic order parameters 

In summary, we have reported the first direct evidence for a variation of the density of 
states of nodal quasiparticles due to a field-induced Doppler energy shift through a four- 
fold field-angle oscillation in the heat capacity. The angular variation, together with the 
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Vll variation of the heat capacity at 2.5 K, strongly suggests that the gap function is 
extremely anisotropic or, more likely, that YNi 2 B 2 C is an unconventional superconductor. 
The magnetic field-angle, field-intensity and temperature dependence of the heat capac- 
ity are in quantitative agreement with a model in which the fourfold pattern arises from 
Doppler-enhanced, fully 3D nodal quasiparticles with momenta in < 100 > directions. 
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FIG. 1: Field dependence of the specific heat at 2.5 K along (100) direction. The open circles are 
with increasing magnetic field and the cross are with decreasing field. The dashed line describes 
C + b(H -# ) 1/2 with C = 23, b = 10. H is essentially due to Meissner effect and 0.09 T is used 
in the fit. The inset is the temperature dependence of the heat capacity around superconducting 
temperature at zero field. 



9 



34.8 



o 




angle 

FIG. 2: Field directional dependence of the heat capacity at 2 K in 1 T. The field angle a is 
measured with respect to the a-axis. The top panel shows the total heat capacity (open circles) 
and the 2-fold component, Ciipt) relative to baseline of 34.1 (dashed line). The bottom panel shows 
the same data after subtracting the background, C±{a) = C tota i(a) — Co — (72(a). The solid line 
describes a fit with a cusped function, C±(a) = 04(1 + T| sin2a|) with T = 0.05. 
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FIG. 3: The 3D surface plots describe magnetic field-angle (x-axis) and field-intensity (y-axis) 
dependence of the in-plane specific heat ratio at 2 K: C(h,a)/C(h,0). Fig. (a) is a numerical 
calculation of the specific heat due to nodal quasiparticles in 3D system and Fig. (b) is experimental 
data taken at 2 K in vortex phase. The experimental data were FFT low-pass filtered with a cut-off 
frequency of 0.03Hz with angles being considered as time. Both plots (a) and (b) were constructed 
with 0.6, 1, 1.5, 2 T data sets. 
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FIG. 4: The oscillation amplitude r at 6 = 90° is shown against the reduced field, H/H C 2 at 2 K 
where H C 2 is 6.3 T. The dashed lines in the main graph and in the inset are numerical estimations 
of the angular contrast by nqp's in 3D with the Fermi velocity vf = 1.3 x 10 7 cm/sec. Inset: The 
oscillation amplitude was shown against temperature in 1 T. 
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